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SUSY-QCD corrections in the squark-gluino sector 
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A status report is given of the calculations of next-to-leading-order (N = f ) supersymmetric QCD corrections 
to the production of squarks and gluinos in pp/pp collisions. The implementation of these SUSY-QCD corrections 
leads to more stable theoretical predictions and to a substantial increase of the production cross-sections. In 
addition we give a discussion of the use of the MS scheme for renormalizing the coupling constants in the QCD 
sector of (iV = 1) supersymmetric theories. 
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1. INTRODUCTION 

The colored squarks ((1l,Qr) and gluinos (g), 
the supersymmetric partners of the quarks (q) 
and gluons (g), can be searched for most effi- 
ciently at high-energy hadron colliders. As R- 
parity is conserved in the QCD sector of (N = 1) 
supersymmetric theories, these particles are al- 
ways produced in pairs. At the moment they can 
be searched for at the Fermilab Tevatron, a pp 
collider with a centre-of-mass energy of 1.8 TeV. 
In the future the CERN Large Hadron Collider 
(LHC) , a pp collider with an envisaged centre-of- 
mass energy of 14 TeV, will allow to cover mass 
values up to 1-2 TeV. 

At present the main search strategy employed 
at the Tevatron involves the search for jets and 
missing transverse energy jl],|| . R-parity conserv- 
ing decay cascades of the squarks and gluinos re- 
sult in Standard-Model hadrons, i.e., jets, and 
LSP's (x?), responsible for the missing energy. 
An example of a complete process, involving pro- 
duction and subsequent decay, is given by 
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PP^ q~L + q R ^ X + q + XiQ 
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with q' the isospin partner of q. The present sta- 
tus of the squark and gluino searches, based on 
this strategy, is depicted in Fig. |]. As the search 
has been negative up to now, only lower bounds 
on the squark and gluino masses are given. The 
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Figure 1. Experimental lower mass bounds for 
squarks (q^t) and gluinos from the Tevatron [Q. 



present lower mass bounds, derived from Fig. [|, 
are (q ^ t) 

m~ g > 175 GeV, 

rriq > 175 GeV for m~ g < 300 GeV. 

So far this analysis has been based on the 
lowest-order (LO) production cross-sections. For 
obtaining adequate theoretical predictions the LO 
cross-sections are in general not sufficient. The 
most important arguments in favor of an analysis 
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that takes into account the next-to-leading-order 
(NLO) SUSY-QCD corrections are: 

• The LO cross-sections have a strong depen- 
dence on the a priori unknown renormal- 
ization scale Qr. Consequently, the theo- 
retical predictions have in general an uncer- 
tainty that is almost as large as the cross- 
section itself. By implementing the NLO 
corrections a substantial reduction of the 
scale dependence is expected. 

• From experience with similar processes 
(e.g. hadroproduction of top quarks), the 
NLO QCD corrections are expected to be 
of the order of +30%. 

• An enhancement of the cross-section would 
lead to a higher value for the lower mass 
bounds for squarks and gluinos. 

• In case of discovery of squarks and gluinos, 
a precise knowledge of the total cross- 
sections is required for the determination of 
the masses of the particles. In contrast to 
the production of e.g. top-quark or Z-boson 
pairs, it seems unlikely that the masses of 
squarks and gluinos can be determined by 
means of reconstruction (in view of the in- 
visible LSP's). Probably they can only be 
inferred from the experimental production 
rates and the theoretical prediction for the 
total cross-section. 

Here we report on the status of the calculations 
of NLO SUSY-QCD corrections to the production 
of (on-shell) squarks and gluinos in pp/pp colli- 
sions, based on the studies presented in j|. For a 
more comprehensive report on this topic we refer 

to 1. 

2. TECHNICAL SET-UP 

We consider the following hadronic production 
processes (generically depicted in Fig. ||): 

pp/pp -> qq, qq, gg, qg (q ^ i), (1) 

where the chiralities and flavors of the squarks 
(e.g. ul, dn) as well as the charge-conjugate fi- 
nal states (e.g. qq) are implicitly summed over. 




Figure 2. Generic Feynman diagram for the pro- 
duction of squarks and gluinos in pp/pp collisions. 



In analogy to the experimental analysis, we ex- 
clude top-squarks from the final state and take 
all produced squarks to be mass degenerate. A 
study of the production of pairs of top-squarks 
is in progress, including the (model-dependent) 
mixing effects in the squark sector. At the par- 
tonic level (right-hand side of Fig. |J) many differ- 
ent subprocesses contribute at LO and NLO, cor- 
responding to different flavors/chiralities of the 
squarks and different initial-state partons. The 
initial-state partons are made up of the massless 
gluons and the five light quark flavors (n/ = 5), 
which are considered to be massless as well. Note 
that not all initial states are possible for a given 
final state. At LO, for instance, the production 
of squark-antisquark final states requires quark- 
antiquark or gluon-gluon initial states, whereas 
the squark-gluino final states are only possible in 
quark-gluon reactions. 

The NLO SUSY-QCD corrections comprise the 
virtual corrections (consisting of self-energy cor- 
rections, vertex corrections, and box diagrams), 
real-gluon radiation (with an additional gluon at- 
tached to the LO diagrams), and the radiation of 
a massless quark (opening additional initial-state 
channels: e.g. gq — > qqq). 

For the particles inside the loops we use the 
complete supersymmetric QCD spectrum, i.e., 
gluons, gluinos, all quarks, and all squarks. We 
have excluded the top-squarks from the final 
states. In order to have a consistent NLO calcula- 
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Figure 3. Feynman diagrams for the virtual NLO SUSY-QCD corrections to the quark-squark-gluino 
vertex (Yukawa coupling) and the quark-quark-gluon vertex (gauge coupling) . 



tion, however, we have to take these top-squarks 
into account inside loops. For the sake of simplic- 
ity we take them to be mass-degenerate with the 
other squarks. For the top quark we use the mass 
m t = 175 GeV. Consequently the final results will 
depend on two free parameters: the squark mass 
uiq and the gluino mass rrig. 

For the internal gluon propagators we use the 
Feynman gauge. As a result, Faddeev-Popov 
ghost contributions have to be taken into account 
in the gluon self-energy and the three-gluon ver- 
tex corrections. For the external gluon lines only 
the transverse polarizations are needed. 

The divergences appearing in the NLO cor- 
rections are regularized by performing the cal- 
culations in n = 4 — 2e dimensions. These di- 
vergences consist of ultraviolet (UV), infra-red 
(IR), and collinear divergences, and show up as 
poles of the form (i — 1,2). For the treat- 
ment of the 75 Dirac matrix, entering through the 
quark-squark-gluino Yukawa couplings, we use 
the 'naive' scheme. This involves a 75 that anti- 
commutes with the other gamma matrices. This 
is a legitimate scheme at the one-loop level for 
anomaly-free theories. 

The UV divergences can be removed by renor- 
malizing the coupling constants and the masses 
of the massive particles. In the case of the mass 
renormalization we have opted for an on-shcll 
scheme with real subtraction point, involving the 
subtraction of the real part of the on-shell self- 



energies at the real- valued pole masses. For the 
renormalization of the QCD coupling constant 
one usually resorts to the modified Minimal Sub- 
traction (MS ) scheme. The MS scheme involves 
n-dimensional regularization, i.e., treating fields, 
phase space, and loop momenta in n dimensions. 
The UV 1/e poles are subtracted, together with 
specific transcendental constants, at an a pri- 
ori arbitrary subtraction point (renormalization 
scale) Qr . In supersymmetric theories, however, 
a complication occurs. In n ^ 4 dimensions the 
MS scheme introduces a mismatch between the 
number of gluon (n — 2) and gluino (2) degrees 
of freedom. As these O(e) mismatches will re- 
sult in finite contributions, the MS scheme vio- 
lates supersymmetry explicitly. In particular, the 
qqg Yukawa coupling g g , which by supersymme- 
try should coincide with the qqg gauge coupling 
g s , deviates from g s by a finite amount at the one- 
loop order. Requiring the physical amplitudes to 
be independent of the renormalization scheme, a 
shift between the bare Yukawa and gauge cou- 
plings must be introduced in the MS scheme, 



9s = g 3 















(§"- 




= 9s 















(2) 



which effectively subtracts the contributions of 
the false, non-supersymmetric degrees of freedom 
(also called e scalars) ||. Here we used a s = 
<rf/(47r), N = 3, and C F = (N 2 - 1)/(2N) = 4/3. 

The need for introducing a finite shift is best 
demonstrated for the effective (one-loop cor- 
rected) Yukawa coupling (see Fig. @) , which must 
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be equal to the effective gauge coupling in an 
exact supersymmetric world with massless glu- 
ons/gluinos and equal-mass quarks/squarks. For 
the sake of simplicity we define the one-loop cor- 
rected effective couplings Tf(Q 2 ) = g s [1+8(Q 2 )] 
and f f(Q 2 ) = g s [1 + 5(Q 2 )} in the limit of 
on-shell quarks/squarks and almost on-shell glu- 
ons/gluinos, with virtuality Q 2 <C rrv- = m 2 . In 
this limit the couplings do not contain gauge- 
dependent terms. In the M S scheme we find after 
charge renormalization || 



MS: tfm 



rf(Q 2 



+ 



4tt V2 



C F --N 



(3) 



The difference between the two effective couplings 
coincides with the shift Taking into ac- 

count this finite shift of the bare couplings in 
the MS scheme, both effective couplings become 
identical at the one-loop level. In this way su- 
persymmetry is preserved and the MS scheme 
becomes a workable one. 

An alternative renormalization scheme is the 
modified Dimensional Reduction (DR ) scheme. 
This scheme consists in treating the fields in 4 di- 
mensions and the phase space and loop momenta 
in n dimensions. As such no mismatch is intro- 
duced and supersymmetry is preserved. At the 
level of the above-defined effective couplings this 
is reflected in the equality || 



DR 



rf(Q 2 



(4) 



As a result, both couplings are identical order by 
order. It should be noted that the transition from 
the effective gauge coupling in MS to the one in 
DR involves a well-known finite renormalization 

(a s A)/(24vr) =a a /{8n). 

In the following we use the MS renormalization 
scheme in combination with the finite shift of the 
Yukawa coupling. In this way supersymmetry is 
preserved on the one hand, while on the other 
hand the definition of the strong gauge coupling 
corresponds to the usual Standard-Model mea- 
surements. In addition to the poles also some log- 
arithms are subtracted in order to decouple the 
massive particles (top quarks, squarks, gluinos) 
from the running of a s (Q 2 R ). The Q 2 R evolution 
of the strong coupling is in this decoupling renor- 
malization scheme completely determined by the 
light-particle spectrum (gluons and rif = 5 mass- 



less quarks): 
91og(Q«) 



11 *r 2 

— N n r 

3 3 1 



(5) 



The above described methods to renormalizc 
the UV divergences result in cross-sections that 
are UV finite. Nevertheless there are still diver- 
gences left. The IR ones will cancel in the sum 
of virtual corrections and soft-gluon radiation. In 
order to separate soft from hard radiation a cut- 
off A is introduced in the invariant mass corre- 
sponding to the radiated gluon and one of the 
produced massive particles. If soft and hard con- 
tributions are added up, any A dependence disap- 
pears from the cross-sections for A — > 0. The re- 
maining collinear singularities, finally, can be ab- 
sorbed into the renormalization of the parton den- 
sities, carried out in the MS mass-factorization 
scheme. This introduces yet another a priori un- 
known scale, the factorization scale Qf- 

In the context of mass factorization an inter- 
esting observation was made during the compar- 
ison of the MS and DR results. Even after cor- 
recting for e-scalar contributions, leading to fi- 
nite differences between the splitting functions in 
MS and DR (t]], discrepancies of £)(m|,m|) per- 
sisted. The exact source of these discrepancies is 
under investigation. This phenomenon is, how- 
ever, not an artefact of supersymmetry, since it is 
also observed in the (pure QCD) process of top- 
quark production Q . 

The last singularity we have to deal with is of 
kincmatical nature. If the gluinos are lighter than 
the squarks, on-shell squarks can decay into mass- 
less quarks plus gluinos (q — > qg). This means 
that inside the phase space of massless-quark- 
radiation processes like gq — » ggq (see Fig. ||) 
explicit particle poles can emerge. These parti- 
cle poles, however, correspond to on-shell (LO) 
production of the intermediate on-shell squark 
state, with subsequent (LO) decay into a mass- 
less quark plus gluino. In our narrow- width ap- 
proach, consisting in neglecting the finite decay 
widths of the squarks/gluinos whenever possible, 
these situations are already accounted for by the 
LO cross-sections (e.g. squark-gluino production 
in Fig. |i|) . In order to avoid double counting these 
kincmatical situations have to be subtracted from 
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Figure 4. Example of a NLO Feynman diagram 
that can give rise to an on-shcll intermediate 
squark state. 



the NLO corrections. Similar subtractions have 
to be performed if the squarks are lighter than 
the gluinos. For more details concerning the ex- 
act subtraction procedure we refer to Q . 

3. RESULTS 

At the partonic level at least three sources of 
large corrections can be identified. Indicating the 
velocity of the produced heavy particles in their 
centre-of-mass system by (3, two of these sources 
reside in the region near threshold (/3 <c 1), from 
which an important part of the contributions to 
the hadronic cross-sections originates. First of all, 
the exchange of (long-range) Coulomb gluons be- 
tween the slowly moving massive particles in the 
final state (see first generic diagram of Fig. ||) 
leads to a singular correction factor ~ na s /(3, 
which compensates the LO phase-space suppres- 
sion factor /3. It should be noted, however, that 
the finite lifetimes of the squarks/gluinos reduce 
(screen) this effect considerably. Secondly, as a 
result of the strong energy dependence of the 




Figure 5. Generic diagrams leading to the 
Coulomb singularity, the large threshold loga- 
rithms, and the high-energy plateau, respectively. 



cross-sections near threshold, large "soft" correc- 
tions ~ log\p 2 ) (i = 1,2) are observed in the 
initial-state gluon-radiation contribution (see sec- 
ond generic diagram of Fig. ^J). At high energies, 
finally, the NLO partonic cross-sections behave 
asymptotically as a constant, rather than scaling 
with 1/s like the LO cross-sections. This is caused 
by the presence of almost on-shell IR gluons in 
space-like propagators, as appearing in the con- 
tributions from hard-gluon/quark radiation (see 
third generic diagram of Fig. |J) . All these effects 
can be calculated analytically in NLO, providing 
powerful checks on the partonic results . 

After convolution of the partonic cross-sections 
with the relevant parton densities, the hadronic 
cross-sections are obtained. When discussing 
LO and NLO results, we calculate all quanti- 
ties \a s {Q 2 R ) : the parton densities, and the par- 
tonic cross-sections] in LO and NLO, respectively. 
Bearing this in mind the hadronic results can be 
summarized as follows. 

(i) As is exemplified in Fig. ^, we find that the 
theoretical predictions for the production cross- 
sections are stabilized considerably by taking into 
account the NLO SUSY-QCD corrections. In all 
processes, for both the Tevatron and the LHC, it 
is observed that the dependence on the renormal- 
ization/factorization scale Q (Q — Qr — Qf) is 
quite steep and monotonic in LO, whereas the Q 
dependence is reduced by roughly a factor of 2.5- 
4 in NLO for reasonable variations of the scale. 
Even a broad maximum develops for scales that 
are roughly a factor of 3-4 smaller than the av- 
erage mass of the final-state particles. The vari- 
ation of the cross-sections as a result of different 
NLO parametrizations of the parton densities is 
< 10% at the Tevatron and < 13% at the LHC, 
where the gluon densities are more important. It 
should, however, be noted that considering such 
a sample of different parton densities could lead 
to an underestimation of the actual theoretical 
uncertainties. 

(ii) From now on we use GRV94 parton den- 
sities and conservatively take as default scale Q 
the average mass of the produced particles. The 
K- factors, K = <jnlo / vlo , depend strongly on 
the process. This is exemplified in Fig. [?] for the 
Tevatron. For both collider types the NLO cor- 
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Figure 6. The scale and parton-density depen- 
dence of the LO and NLO cross-sections for 
gluino-pair production at the Tevatron. Par- 
ton densities: GRV94 (solid), CTEQ3 (dashed), 
and MRS A' (dotted). Mass parameters: rriq = 
280 GeV, rrig = 200 GeV, and m t = 175 GeV. 



rections are positive and large (up to +90%) for 
the dominant production cross-sections, involving 
at least one gluino in the final state. The corre- 
sponding K-factors also exhibit a sizeable depen- 
dence on the squark and gluino masses (especially 
the rrig dependence of Kg g is large). The NLO 
corrections for squark final states are moderate 
( < + 30%) . In view of the direct link between 
the total cross-sections and the experimental de- 
termination of the squark / gluino masses in case of 
discovery, the inclusion of the NLO SUSY-QCD 
corrections is indispensable. 

(Hi) Comparison of the NLO cross-sections 
with the cross-sections used in the Tevatron anal- 
ysis (LO, EHLQ parton densities, and a scale 
Q that equals the partonic centre-of-mass energy 
fl|,§|) reveals that the NLO corrections raise the 
lower mass bounds for squarks and gluinos by 
+10 GeV to +30 GeV. 

(iv) Apart from the total cross-sections, also 
distributions with respect to the rapidity (y) and 
transverse momentum (p t ) of one of the outgoing 
massive particles can be studied. The K-factors 
for these distributions are independent of y for 
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Figure 7. The K-factors for the Tevatron. Parton 
densities: GRV94 with Q = average mass. Mass 
parameters: m g = 200 GeV and m t — 175 GeV. 



all practical purposes and hardly depend on p t , 
except for large values of pt where the NLO cor- 
rections make the distributions somewhat softer. 
Consequently, multiplication of the LO distribu- 
tions with the above-defined K-factors for the to- 
tal cross-sections approximates the full NLO re- 
sults quite well. 
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